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A simple discussion on the backreaction of inhomogeneities in cosmology, focusing on the 
possibility that it could explain the present acceleration and solve the coincidence problem. 

1 Introduction 

The coincidence problem. Perhaps the most surprising observation in modern cosmology 
is that the expansion of the universe has apparently started to accelerate in the recent past, 
at a redshift of probably less than one. This phenomenon has usually been interpreted in the 
context of a homogeneous and isotropic model of the universe. The acceleration can then be 
accommodated only by modifying gravity or by introducing an energy component with negative 
pressure (or energy density), 'dark energy', the magnitude of which is fit to the observations. 

These models suffer from the coincidence problem: they do not explain why the energy 
densities of matter and dark energy have become comparable only recently. If there is a dy- 
namical explanation for the acceleration having started in the near past -instead of it being a 
coincidence- then presumably it is related to the dynamics we see in the universe. The most 
significant change in the universe at small redshifts is the formation of large scale structure, 
so it seems a natural possibility that the observed deviation from the natural prediction of 
the homogeneous and isotropic model could be related to the growth of inhomogeneities in the 
universe. 

The fitting problem. The reasoning behind using a homogeneous and isotropic model is 
that the universe appears to be homogeneous and isotropic when averaged over sufficiently large 
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scales. In the usual approach one first takes the average of the metric and the energy-momentum 
tensor, and then plugs these smooth quantities into the Einstein equation. However, physically 
one should first plug the inhomogeneous quantities into the Einstein equation and then take the 
average. Because the Einstein equation is non-linear, these two procedures are not equivalent: 

{G^A9a(i))^G,,{{g^p)) , (1) 

where G^y is the Einstein tensor, gap is the metric and () stands for averaging. The content of 
is that the average behaviour of an inhomogeneous spacetime is not the same as the behaviour of 
the corresponding smooth spacetime. Here "corresponding" means that the smooth and average 
quantities have the same initial conditions. In other words, the average properties of an inhomo- 
geneous spacetime (energy density, expansion rate, . . . ) do not satisfy the Einstein equation. 

This is the fitting problem discussed by George Ellis in 1983121. Simply put, how does one 
find the average model which best fits the real inhomogeneous universe? The difference between 
the behaviour of average and smooth quantities is also known as backreaction. 

The equations satisfied by the average quantities in a general inhomogeneous spacetime have 
been obtained!^. However, these equations are not closed, which simply means that different 
spacetimes sharing the same initial average values evolve differently even as far as their averages 
are concerned - not a surprising result. We have taken (followingl^ a more modest approach, 
in which one assumes a smooth background and studies the effect of perturbations. 

2 The backreaction calculation 

The expansion rate. We consider a homogeneous and isotropic universe that is spatially 
flat and filled with a pressureless fluid. We assume that the fluid has scale-invariant adiabatic 
perturbations (with zero mean) and an amplitude given by the measurements of the CMB. In 
first order perturbation theory, the metric in the longitudinal gauge reads 

ds'^ = -{l + 2^{x))dt'^ + {1-2^ {x))a{tfdx'^ , (2) 

where a = {t/t^f'^^ and <I> is constant in time. 

We want to find out the effect of the perturbations on the expansion rate measured by 
a comoving observer. Since the spacetime is inhomogeneous, the expansion rate is different 
for observers in different points, and to make the comparison to a homogeneous and isotropic 
universe, we should take an average. The expansion rate is a covariantly defined scalar quantity, 
so this operation is we ll de fined, but the result depends on which hypersurface one chooses 
to take the average onEEl. We define the average expansion rate to be the average of the 
expansion rates measured by comoving observers (weighed by the volume element) whose clocks 
show the same proper time r (measured from the big bang). In other words, we average on the 
hypersurface of constant proper time r of comoving observers. 

The covariant definition of the local expansion rate is ^ = u^.^, where u^^ is the four-velocity 
of the matter fluid. From the metric 0, we find and from that we obtain the expansion 
rate 6 and the proper time r. Then we simply express 9 in terms of r and take the average. 
Note that this calculation involves no dynamics: we just take a given metric and calculate the 
observable of interest for that spacetime. 

Naively, one might expect the average expansion rate to be 

{6) ~ 3i/,(l + ai($)+ 02 + 03(^2)) , (3) 

where H^- = 2/(3t) is the background expansion rate in terms of the proper time, and Oj are 
some constants of order one. Here and in what follows we expand only to second order in $ (note 



that we use results from first-order perturbation theory, so that this is not a consistent second 
order calculation). Since ^ has zero mean, the linear term and its square vanish, and since the 
amplitude of $ is 10~^, the quadratic term is of the order of 10"^*^ ~ completely negligible. 
When one actually calculates the expansion rate, the result is not ©, but instead 

= 3Hr {l + Xittr + X2al) , (4) 

where Pi are given constants of order one, a,- = (t/tq)^/^, and Aj are constants. 

The plain powers of <I> that one would expect do not appear, and instead the backreaction 
is given by gradients of ^. The structure of the backreaction terms in is easy to understand: 
each gradient has to be accompanied by 1/a for reasons of covariance and by l/-ff for reasons 
of dimensionality. Since is a scalar, gradients have to appear in pairs, and since the Einstein 
equation is second order, there are at most two gradients for each $. Since l/{aH)'^ oc a, the 
backreaction grows relative to the background, and thus behaves like dark energy. 
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Backreaction as dark energy. For a homogeneous and isotropic spatially flat universe, the 
expansion law is (in this case 9 = 3H) 

1 

^l 

where Mpi is the Planck mass and p is the energy density, which for matter behaves like p cx a~^. 
In contrast to ©, the average of the inhomogeneous expansion law is, from 

oc a-3 + 2Aia-2 + (A? + 2A2)a7^ + 2AiA2 + Aia^ , (6) 

where in the second equality we have taken into account oc a~^. 

An observer in a perturbed universe trying to fit the smooth model © to the observed 
expansion rate © would conclude that there is a mysterious energy component which is nowhere 
to be seen but which affects the expansion rate. Writing p = Pm + Pde where pm is the energy 
density of matter, and pde is the apparent energy density of this 'dark energy', one finds from 
the usual relation pd^, oc a^'^^^+'^) the equations of state w = —1/3, —2/3, —1 and —4/3. 

The equation of state is negative because the backreaction grows relative to the background 
for which w = 0. Note that there is nothing unnatural about having an equation of state which is 
more negative than —1 (as marginally indicated by observations ). It does not imply violation 
of the weak energy condition, since the 'dark energy' is only a parametrisation of our ignorance 
of the real expansion law, and does not correspond to an actual energy component. 

Backreaction has a negative equation of state that could produce acceleration, so the next 
question is whether the effect is significant. In other words: what are the values of Ai and A2? 
For perturbations in the linear regime, we have 

Ai = [3i{di^d^^)/Hl ~ 10-5 , A2 = (i2{^^{V^m^))/H^ - {5l) ~ 1 , (7) 

where Hq is the Hubble parameter today and 5q is the density perturbation today. The constant 
Ai is easy to evaluate and is negligible, while A2 is more involved. We have used the relation 
= -3{aHf'5/2, valid in the linear regime of perturbations, which we have taken to end at 
(5^) = 1. The quantity to be averaged can then (by definition) be of order one point by point. 
However, since it is a total derivative, the periodicity implied by the standard decomposition in 
terms of Fourier modes puts it artificially to zero. The contribution of similar terms at order <I>^ 



could be non-vanishing, since the square of a total derivative is not a total derivative. If we 
would have taken the average of (P' instead of averaging Q and squaring it, we would have found 

((^/3)2) ^ /^2(^i + ±_l_(vW$))=/72(^l + i(52)) , (8) 

where all subdominant terms have been dropped. For comparison to the observations it is more 
correct to take the average first, but Q shows that backreaction from linear perturbations can 
have a sizeable effect; taking the cut-off to be at ((5^)=1, the correction is of the order of 10%. 

The reason that the effect can be large even though the metric perturbation $ is small is 
that the physics is not given just by the metric but by its first and second derivatives. The 
derivatives of the metric are dimensional quantities, so they are small or large only with respect 
to some scale. In the case of spatial derivatives in cosmology, the comparison scale turns out to 
be a/f, the size of the horizon, compared to which the spatial gradients are large. 

3 Conclusion 

In summary, the backreaction of linear perturbations on the expansion rate is non-zero, boosted 
by powers of k"^ /{aH)"^ from the naive expectation of powers of and naturally involves a 
negative equation of state. For periodic boundary conditions, it is also numerically negligible at 
second order in 

In evaluating the backreaction of linear perturbations, we had to introduce a cut-off to obtain 
a finite result. A more realistic treatment would take non- linear perturbations and structure 
formation explicitly into account; note that the presence of 5 relates backreaction directly to the 
growth of structure. This calculation has not been done, but one would expect the growth of 
the gravitational potential to increase the magnitude of backreaction and make its equation of 
state more negative. The backreaction of perturbations breaking away from the linear regime 
could thus possibly give 'dark energy' at the right time to solve the coincidence problem. 

Whether or not backreaction turns out to be quantitatively important, the impact of in- 
homogeneities on the expansion rate should be properly evaluated to solve the fitting problem 
and make sure that we are using the right theoretical equations to explain increasingly precise 
cosmological observations. 
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